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Abstract 

We propose a unitary matrix Chern-Simons model representing fractional 
quantum Hall fluids of finite extent on the cylinder. A mapping between 
the states of the two systems is established. Standard properties of Laugh- 
lin theory, such as the quantization of the inverse filling fraction and of the 
quasiparticle number, are reproduced by the quantum mechanics of the ma- 
trix model. We also point out that this system is holographically described 
in terms of the one-dimensional Sutherland integrable particle system. 
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1 Introduction 



Chern-Simons actions for gauge fields, since their introduction to physics |l|], have 
found numerous apphcations in elementary particle and condensed matter situations. 
In the last of these developments, Chern-Simons theory on the noncommutative 
plane has been proposed by Susskind as an effective description of the fractional 
quantum Hall fluid . Specifically, the ground state of this theory can be interpreted 
as the Laughlin state for an infinite number of electrons 0. The filling fraction 
corresponds to the inverse coefficient of the Chern-Simons action. 

The basic idea behind this identification is that the noncommutative Chern- 
Simons action describes a (noncommutative) magnetic membrane, which, in turn, 
is equivalent to a magnetic fluid. The connection between (commutative) mem- 
branes and fluids is old § and has recently been extended to supersymmetric 
fluids with spin 0. The new ingredient in the quantum Hall connection is the pro- 
posal by Susskind that an essentially noncommutative fluid is appropriate in order 
to incorporate the discrete nature of electrons. This noncommutativity exists al- 
ready at the classical level, and is distinct from the (quantum) noncommutativity 
of the coordinates of particles in the lowest Landau level ^ ^ and its extension to 
magnetohydrodynamics . 

The above Chern-Simons theory can describe only an infinite number of electrons 
on an infinite plane. In a previous paper we proposed a regularized version of the 
noncommutative theory on the plane in the form of a Chern-Simons matrix model 
with boundary terms . This model describes a system of finitely many electrons 
(a quantum Hall 'droplet') and was shown to reproduce all the relevant physics of 
the finite Laughlin states, such as boundary excitations, quantization of the filling 
fraction and quantization of the charge of quasiparticles (fractional holes). We 
further pointed out that the matrix model, and thus also the quantum Hall system, 
is equivalent to the Calogero model [^-[[15], a one-dimensional system of particles 



whose connection to fractional statistics [|I6|, anyons [|18|, |T9|, ^ and the 



quantum Hall system 23| has long been established. An explicit (although non- 
unitary) mapping between the states of the matrix model and Laughlin states was 



presented in [24 



It is of interest to extend the correspondence between the noncommutative 
Chern-Simons matrix model and the quantum Hall system for spaces of different 
topologies which have compact dimensions. This is of both theoretical and practical 
significance, since compact spaces provide a natural regularization and have been 
used in alternative approaches to Laughlin states. 

In this paper we shall present such a generalization for the case of a space with 
cylindrical topology, in the form of a unitary matrix model. We shall identify and 
analyze its classical and quantum states and establish a mapping to Laughlin states. 
Similarly to the planar case, we shall demonstrate that this model is equivalent to 
a periodic version of the Calogero model known as the Sutherland model. Finally 
we shall conclude with some directions for future research. 
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2 Chern-Simons theory on the noncommutative 
plane and quantum Hall states 



Before analyzing the problem on the cylinder, we will review the basic features of 
Chern-Simons (CS) theory on a noncommutative plane and a commutative time 
and its connection with quantum Hall states, as proposed by Susskind 0. We will 
also briefly review the corresponding Chern-Simons matrix model describing flnite 
quantum Hall droplets |I0| . 



2.1 Noncommutative Chern-Simons theory from magnetic 
fluids 

The system to be described consists of an incompressible fluid of oo spinless 

electrons on the plane in an external constant magnetic fleld B (we take their charge 
e = 1). Their coordinates are represented by two (inflnite) hermitian matrices X^, 
a = 1,2, that is, by two operators on an inflnite Hilbert space. The average electron 
density is po = ^/^ttO. The action is the analog of the gauge action of particles in a 
magnetic fleld: 

S = Jdt^Tr {eabiXa + t[Ao, Xa])X, + 2dAo} (1) 

with Tr representing (matrix) trace over the Hilbert space and [., .] representing 
matrix commutators. Similar actions flrst appeared in matrix Chern-Simons the- 
ory as a possible approach to the fundamental formulation of M-theory [^. The 
speciflc action above has the form of a noncommutative CS action in the operator 
formulation [^]. Gauge transformations are conjugations of Xa by arbitrary time- 



dependent unitary operators which are compact enough to leave traces invariant. 
These have nontrivial topology and lead to level quantization (see p8| for an 
analysis of this class of transformations). In the quantum Hall context they take 
the meaning of reshuffling the labels of the electrons, a generalization of particle 
permutation operators. Equivalently, the Xa can be considered as coordinates of 
a two-dimensional fuzzy membrane, 2n6 playing the role of an area quantum and 
gauge transformations realizing area preserving diffeomorphisms. 

The time component of the gauge fleld ensures gauge invariance, its equation of 
motion imposing the Gauss law constraint 

-tB[Xi,X2] = B9 = ^ (2) 

with u = 2ttpq/B being the fllling fraction. The canonical conjugate of Xi is 
P2 = BX2, so the operator in the left-hand side of (0) is the generator of gauge 
transformations on Xa- Since gauge transformations are interpreted as reshufflings 
of particles, (^ above has the interpretation of endowing the particles with quantum 
statistics of order 1/z/. 
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We will assume that Xi,X2 provide an irreducible representation of the Gauss 
law (|]), else we would be describing multiple layers of quantum Hall fluids. This 
representation is essentially unique, modulo gauge transformations, so there is a 
unique state in this theory (the vacuum). Deviations from the vacuum state can be 
achieved by introducing sources in the action p|. A localized source at the origin 
has a density of the form p = pq — q6^{x) in the continuous (commutative) case, 
representing a point source of particle number —q, that is, a hole of charge q for 
g > 0. The noncommutative analog of such a density is 

[X,,X2]=i9{l + q\0){0\) (3) 

where \n), n = 0, 1, . . . is an oscillator basis for the (matrix) Hilbert space, |0) 
representing a state of minimal spread at the origin. In the membrane picture the 
right-hand side of (^ corresponds to area and imphes that the area quantum at the 
origin has been increased to 271$ {1 + q), therefore piercing a hole of area A = 27i6q 
and creating a particle deficit q = poA. We shall call this a quasihole state. For 
g > a solution of (0) is 

oo 

X^ + iX2 = V2eJ2y/WTq\n-l){n\ (4) 

n=l 

The above assumes that |0) is really a state at the origin, meaning (Xi +2X2)10) = 0. 
Without this residual condition (|) has many more solutions, since |0) is a gauge- 
dependent state that can be reshuffled around. For instance, a class of solutions 
is 

°° , 

Xi + iX2 = P V29 ^ + qd{n - Uo) \n - 1) {n\ P (5) 

n=l 

with ■d{s) the usual step function {^{s) = 1 if s > 0, else 'd{s) = 0) and P the 
permutator of |0) and \no) 

p = i-(|o)-K))((o|-K|) (6) 

This represents an annular hole of charge q at distance ~ ^/On^ from the origin. 

For quasiparticles (g < 0), as long as — g < 1 we have a similar equation and 
solution. For — g > 1, however, clearly equations (^ and (^ cannot hold since the 
area quantum cannot be diminished below zero. The correct equation is, instead, 

[Xi,X2] = ze(i- J2 - ^\k){k\] (7) 

\ n=0 / 

where k and e are the integer and fractional part of the quasiparticle charge — g. 
The solution of (|^) is 

fc— 1 00 

Xi + iX2 = ^n\n) {n\ + V2e ^ Vn - k - e\n - l){n\ (8) 

n=0 n=fc+l 
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where again we assumed that (Xi + iX2)\k) = 0, so \k) now represents the state at 
the origin. In the membrane picture, k quanta of the membrane have 'peeled' and 
occupy positions Zn = Xn + Wn on the plane, while the rest of the membrane has a 
deficit of area at the origin equal to 27r^e, leading to a charge surplus e. The quanta 
are electrons that sit on top of the continuous charge distribution. If we want all 
charge density to be concentrated at the origin, we must choose all Zn = 0, which 
means that Xi + 1X2 annihilates all \n) for n = 0, 1 . . . fc. 

More general quasihole (particle) states, with the holes (or fractional part of the 
particles) positioned at arbitrary points on the plane can easily be constructed, but 
we shall not do so here. We point out that the above particle states for integer q 
are identical to flux solitons of noncommutative gauge theory pOl 



31, 



2.2 Finite Chern-Simons matrix model 

For a finite number of electrons N we take the coordinate Xa to be finite hermitian 
N X N matrices. The action (|l|), however, and the Gauss law (||) are inconsistent 
for finite matrices, and a modified action must be written. We take [IC] 

S = J dt^Ti {eatiXa + i[Ao,Xa])Xh + 26 Ao - ooX^,} + ^^{i^ - Ao^) (9) 
is a complex A^-vector that transforms in the fundamental of the gauge group 

u{Ny. 

Xa ^ UX^U-' , ^^U^ (10) 

while the term proportional to u serves as a spatial regulator providing a harmonic 
potential that keeps the electrons near the origin. 
The Gauss law now reads 

G = -iB [Xi, X2] + - BO = (11) 

Taking the trace of the above equation gives 

^t^ = NB9 (12) 

The equation of motion for in the Aq = gauge is = 0. So we can take it to be 



^ = VNB\v) (13) 

where \v) is a constant vector of unit length. So the traceless part of ([TT]) reads 

[X,,X2]=te{l-N\v){v\) (14) 

This is similar to (^ for the infinite plane case, with an extra projection operator, 
which is the minimal deformation of the planar result (^ that has a vanishing trace. 
\E' clearly acts like a boundary term, absorbing the 'anomaly' of the commutator 



4 



[Xi, X2], very much like the case of a boundary field theory required to absorb the 
anomaly of a bulk CS field theory. 

The classical states of this theory are given by the set of matrices A = Xi + 1X2 
satisfying (0), and can be explicitly found |T^. The ground state is found by 



minimizing the potential V = {Buj/2)TT{Xf + X|) while imposing the constraint 
(p^). We obtain the solution 

N~l 

Xi + «X2 = ^ v^|n- l)(n| , \v) = \N - 1) (15) 

n=0 

This is essentially a quantum harmonic oscillator projected to the lowest N energy 
eigenstates. The radius squared matrix R"^ = Xf + X| has a finite, equidistant 
spectrum. So the above solution represents a circular quantum Hall 'droplet' of 
radius y/2N6 and particle density po = N/{tiR'^) ~ l/(27r^) as in the infinite plane 
case. 

Excitations of the classical ground state can be considered. A class of such 
excitations are perturbations of A = Xi + 1X2 generated by the infinitesimal trans- 
formation 

N-l 

A' = A + ^n{A^r (16) 

n=0 

with e„ infinitesimal complex parameters. The sum is truncated to X — 1 since 
is an X X X matrix and only its first X powers are independent. These map the 
boundary of the droplet to the new boundary (in polar coordinates) 

N 

R'{(j)) = V2N9 + Cne*""^ (17) 

n=-N 

where the coefficients 

Cn = c*_^ = -^^"-1 {n> 0), Co = (18) 

This is an arbitrary area-preserving deformation of the boundary of the droplet, 
truncated to the lowest X Fourier modes. The above states are, therefore, arbi- 
trary area-preserving boundary excitations of the droplet ^ , appropriately 
truncated to refiect the finite, noncommutative nature of the system. 

A second class of excitations are the analogs of quasihole and quasiparticle states. 
States with a quasihole of charge —q at the origin are of the form 

A = V29\j^\N){0\+Y^/irTq\n-l){n\j , q>0 (19) 

representing a circular droplet with a circular hole of area 27r6'g at the origin, that is, 
with a charge deficit q. Note that (|19D stills respects the Gauss constraint (|l^) (with 
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\v) = — 1)) without the exphcit introduction of any external source. The hole 
and the boundary of the droplet together cancel the anomaly of the commutator, 
the outer boundary part absorbing an amount N + q and the inner (hole) boundary 
producing an amount q. 

Fractional quasiparticle states cannot be written in this model, reflecting the fact 
that such states do not belong to the z/ = 1/ B6 Laughlin state. Particle states with 
an integer particle number —q = m and the extra m electrons positioned outside 
the droplet do exist, but we shall not write their explicit form here. 

We conclude by pointing out that boundary excitations, quasiholes and particle 
state can all be continuously deformed to each other, due to the finite number of de- 
grees of freedom of the model. Such transformations become highly nonperturbative 
in the N ^ oo limit. 



3 A model for finite number of electrons on the 
cylinder 

The proposed model works well for electrons on an infinite plane. For a space 
representing a cylinder of radius R we take one of the coordinates, say, X2, to be 
periodic with period 2'irR. Clearly the above model does not take into account 
this periodicity and has to be suitably modified in order to correctly describe the 
physics in the compact dimension. We shall propose here such a model, appropriate 
to describing one compact dimension. 

There are two routes for achieving this. The first is to write a matrix model 



on the covering space and then reduce it [^, |3^, leading to matrices depending on 



an additional continuous parameter dual to the compact dimension, that is, a field 
theory. The second is to represent the compact dimensions with unitary matrices 
P^ . As we shall demonstrate, the two approaches turn out to be equivalent in our 
case. 



3.1 The Chern-Simons unitary matrix model 

We shall begin with the second approach, which is simpler and leads more directly to 
the desired model. The main point is that the coordinate X2 is not single-valued on 
the toroidal space and thus is not a physical observable. An alternative coordinate 
which is single- valued on the cylinder is the exponential e*^^/^. For a noncommuta- 
tive space we define the unitary operator 

U = e'^^l^ (20) 

Together with the hermitian operator Xi = X, they parametrize a noncommutative 
cylinder. The planar noncommutativity relation for Xi,X2 translates into 

UXU-^ = X + ^ (21) 
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To write the Chern-Simons action on such a space we imitate again the magnetic 
action of a particle with coordinates Xj. In the Landau gauge Ai = 0, A2 = —BXi 
the lagrangian reads —BX1X2. Representing X2 as —iRU~^U and including a 
Lagrange multiplier Aq which reproduces (pTF), the full action becomes 



S = J dtBTT[tRU-\U + i[Ao,U])X + eAo} (22) 

It is again expressed in terms of a covariant time derivative DqU of the unitary 
operator U. Note that, in the above, we have adopted the ordering U~^U for 
the operator representing X2. Had we adopted the ordering UU^^ we would have 
ended with an action involving DqU U~^X instead of U~^DqUX . The two actions 
represent identical physics, upon redefining X — > UXU~^ (which, upon use of the 
noncommutativity relation is simply a shift X ^ X + {6 / R)). 

To describe a finite system consisting of N electrons we need to take the coordi- 
nates to be finite N x N matrices. The constraint (|2l|) , however, is not consistent 
for finite matrices, just as in the planar case, and we need to modify the action with 
'boundary' terms that render a consistent form for the constraint. For our purposes 
we take 

S = J dtBTii^tRU-^il + i[Ao,U])X + eAo-'^X'^^ + - Ao^) (23) 



It is similar to the planar Chern-Simons matrix model in [0. The boundary term 
involves \E', a complex X- vector boson that transforms in the fundamental of the 
gauge group U{N). Its role is to absorb the 'anomaly' of the group commutator 
UXU~^ — X, analogous to a boundary field theory required to absorb the anomaly 
of a bulk CS field theory. We also added a spatial regulator term in the form of 
a harmonic oscillator potential in the direction of the cylinder's axis. Since the U- 
direction is compact we need not worry about localizing particles there, and we only 
want to localize them in the infinite X-direction. 

Before analyzing this model further, we present the alternative derivation in 
terms of the covering space reduction and demonstrate that it produces a model 
equivalent to the model proposed above. 



3.2 The Chern-Simons matrix field theory 

An alternative approach to deriving the desired matrix model is to augment the 
dimensionality of the hermitian matrices Xi,X2 of the planar model in [jlOl from 
N to pN and take p 00. p represents the copies of the cylinder on the (planar) 
covering space. To ensure that the state is the same on all copies we must impose the 
condition that fields in different copies are gauge equivalent; that is, the operator 
which shifts copies is a unitary (gauge) transformation. Therefore, there should 
exist some unitary matrix U, representing shifts by one copy, satisfying 

UXiU = Xi , UX2U~^ = X2 + 2ttR 
UAqU-^ = Aq , [/^ = e*"^ (24) 
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This can be explicitly realized by parametrizing the indices /, J of {Xa)ij, {Aq)ij 
and as 

I^i + nN, i^l,...N, n = . . . - 1, 0, 1, . . . (25) 

That is, split Xa in terms of N x N blocks, with the diagonal (n, n) blocks repre- 
senting the electrons on the n-th copy in the covering space and the off-diagonal 
{n, m) blocks representing effective 'interactions' of electrons between the n-th and 
m-th copies on the covering space. Clearly the (n, n) copy must be the same as 
the (0, 0) copy, only shifted by 27rRn in the X2 direction. Further, the interactions 
between the m and n copies must only depend on the distance between the copies 
m — n. Similarly, the Lagrange multiplier {Ao)m.,n must impose constraints on the 
m and n copies that depend only on their distance m — n. Finally ^n, representing 
the boundary of the state in copy n, must be the same for all n, up to an irrelevant 
phase. Overall we have 

{^l)m,n = {^l)m-n , (^2)m,n = {^2)m-n + '^T^Rn^mn 

{Ao)m,n = {Ao)m-n, *n = 6*""* (26) 

The unitary transformation U is simply the shift n ^ n + 1. 

We see that the matrices Xa and Aq are now parametrized in terms of an addi- 
tional integer n. Hermiticity of Xa and Aq in the original indices /, J means 

{Xa)n = (^i)-n , {Ao)n = (4)-n (27) 

We can, therefore, define the Fourier transforms 

Xa{a) = ^(X„)„e-- , Aoia) = ^(A)-^'"" (28) 

n n 

with (7 a variable with periodicity 27r. Xa{(j) and Ao{a) are hermitian cr-dependent 
N X N matrices. Matrix multiplication in the original /, J indices translates into 
matrix multiplication pointwise in a, while /-trace translates into cr-integration and 
matrix trace. It is also useful to define 

^ (a) = ^ ^^e'""" = *5((7 + a) (29) 

n 

We can now write the original matrix model (with a confining harmonic potential 
in the Xi direction) 

ldtBTrl^-X^{X2 + i[Ao,X2]) + eAo-^X^^ + ^\i^-Ao^) (30) 

in terms of the matrices Xi{a),Ao{a). A standard calculation leads to the result 

S = J dtda B Tr |-Xi(X2 + i[Ao, X^] - Rd^A^) - ^X^ + 9Aq^ 

+ j dt *^(ii - Ao((7 = a)*) (31) 
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This is nothing but 1 + 1- dimensional U{N) Yang-Mills theory with a Wilson line 
source at a = a and a uniform background charge. To see this, rename t = ao, 
R^^a = o"i, Xi = —F/uo, X2 = Ai. Then the above action becomes 

S = y rfV ^ Tr (^FFoi - ^F^ + Buj9Ao^ + J dt ^^(z^ - Ao{a = a)^) (32) 

where we defined the field strength 

Foi = doA,~diAo + t[Ao,A,] (33) 

We recognize the Yang-Mills action in the first-order formalism, on a circular space 
of radius R~^. In addition, there is a constant background U{1) charge density Buj6 
and a localized source at ai = a depending on \E'. The latter corresponds to an 
insertion of a Wilson line in the temporal direction carrying the direct sum of all 
symmetric representations of the gauge group U{N). 

To make this last point explicit, consider the temporal direction ctq compact 
and euclidean with period T. This also turns Aq into iA^. An appropriate gauge 
transformation can render AQ{aQ,a) diagonal and independent of ctq. The diagonal 
elements (^o)ii (no sum in j) at ai = a correspond to the eigenvalues of the temporal 
Wilson line element at o"i = a, e*^^ , specifically 

In this gauge the components of \Ef decouple and become N independent bosonic har- 
monic oscillators with frequencies i{Ao)jj. Integrating out \E' produces the partition 
function of these N oscillators which, assuming normal ordering, is 



= 1 + ^ e'^^ + e''^' + E e'^^""'^' + • • • (35) 

3 \ 3 j<k / 

We recognize the terms in the last sum as the trace of the temporal Wilson loop 
in the singlet, fundamental, doubly symmetric etc. representations. The above is, 
then, a (bosonic) oscillator representation of a Wilson loop element, reproducing all 
symmetric representations. A discussion of arbitrary representations in the context 
of the planar CS matrix model can be found in |^ . 



It is a remarkable fact that the above theory can be reduced to a unitary matrix 
model identical to the one derived in the previous section. The details are explained 
in [0 . The basic point is that two-dimensional Yang- Mills theory has no propagat- 
ing modes and the only dynamical degrees of freedom are in the nontrivial holonomy 
of the gauge field Ai around the spatial direction. Defining the U{N) line element 

U[a,b]= Pe'^'^""^' (36) 
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the phase space variables in the gauge Aq = are 



U = U[a, 271 + a] , X = -- /^^^" U[a, a] FU[0, a]-^d(T (37) 



In terms of these the action reduces to 

S = J dtBTi l^iR U-^UX - ^X^j + i^^i! (38) 

The Gauss law for the fields Ai and F involves the background charge 9 and the 
point source at fXi = a. Expressed in terms of the reduced phase space variables 
it becomes the constraint 

X - UXU-' = -^M/vI/t - I (39) 

Inserting the above constraint in the action through a new field Aq (which is now 
only a function of time) we recover the unitary Chern-Simons matrix model (p3|). 
Thus, both routes to compactifying X2 lead to the proposed unitary matrix model. 



4 Classical states of the unitary matrix model 
4.1 General solution 

To get a feeling of the physics of the above model we shall analyze its classical 
structure and states. We can again impose the Aq equation of motion as a Gauss 
constraint and then put Aq = 0. In our case it reads 

G = RB {UXU-^ -X) + -Be = Q (40) 

The trace of the above equation gives as in the planar case 

¥^ = NBO (41) 

In the Aq = gauge \1/ is constant, take it V NB \v) with |t>) a constant unit vector. 
The traceless part of (l40| ) reads 

UXU-'-X = ^il-N\v){v\) (42) 

This is similar to (^) for the infinite cylinder, with an extra projection operator. 
Again, using the residual time-independent U{N) gauge freedom to rotate \v) to the 
form \v) = (0, ... 0, 1) we obtain the form {9/R) diag (1, . . . , 1, 1 — iV) for the above 
constraint, which is the minimal deformation of the result (0) that has a vanishing 
trace. 

The equations of motion for X and U read 

X + [U-'^il, X] = , iRU-^U -ujX = (43) 
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or 

X = i-^{U'^U) = (44) 

This represents free motion on the manifold U{N), as represented by the unitary 
matrix U{t), with X playing the role of matrix momentum. It is solved by 

U{t) = f/oe-'"^°*/^ , X{t) = Xo (45) 



where the constant matrices Xq, Uq satisfy the constraint (^2]). We can find all the 
classical states of this model by diagonalizing Uq = diag {e^'^'^} . By examining the 
diagonal and off-diagonal elements of ( |42D we see that the components of \v) have 
length |t>„p = Choosing their phases as f„ = exp{—i(l)n/2)/\/N, as we can do 

using the residual f/(l)^ gauge invariance, we obtain 

/-n:sm( ^'"^^" ) 

The solution is parametrized by the compact eigenvalues of Uo, (pn, and the 
diagonal elements of X, Xn, that is, by N coordinates on the cylinder. 

4.2 Classical ground state 

The lowest energy state, that is, the state most closely packed around Xi = 0, is 
found by minimizing the potential {Buj/2)X'^ while respecting the constraint (^2|). 
Implementing it with a matrix Lagrange multiplier A we obtain 

BuX + U'^KU - A = , [X, U-'^A.U] = (47) 

The above is solved by 

eBuj(N+l \(N+l \^ . 

Amn = ( nj ( n+ljdmn , Vn = On,l (49) 

The eigenvalues of U are exp{i(j)o + i27m/N) and are evenly distributed on the 
compact dimension of the cylinder. Similarly, the eigenvalues of X are evenly dis- 
tributed along the axis of the cylinder and span a length ~ N6/R. Therefore, the 
above solution represents a tubular quantum Hall droplet around the cylinder with 
an area A = 2ttR ■ N6/R and an average density po = N/A = 1/{2tt6). 

We point out that the average distance between successive electrons in the xi 
direction is 9/R. In the planar case, particles were evenly distributed on the plane 
with a density 1/2t(9 and thus an average distance of order \/d. This is a signal 
that, on the cylinder, quantum Hall states do not have a constant density. In the 
extreme case 6 ^ R^ the electrons will behave more like one- dimensional particles 
with well-defined positions along the length of the cylinder. We also point out the 
existence of 0o in the solution for U, which does not affect the energy. 
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4.3 Equivalence to the Sutherland model 



Just as in the planar case, the matrix model above is equivalent to a one-dimensional 
particle system, the so-called Sutherland model [|12]. This is an integrable system 
of N nonrelativistic particles on the circle with coordinates and momenta {4>n,Pn) 
and hamiltonian 

N —2 

It can be thought of as the Calogero model of particles on the line with inverse- 
square mutual potentials, rendered periodic in space with period 2nR. In terms 
of the parameters of the model, the mass of the particles is B/u and the coupling 
constant of the two-body inverse-square potential is z/~^. We refer the reader to 



T^, |T5|, |5T[] for details of the derivation of the connection between the matrix model 



and the Sutherland model and to [40, 42 for the connection between 2-dimensional 



Yang-Mills and the Sutherland model. Here we simply state the relevant results and 
give their connection to quantum Hall quantities. 

The positions of the Sutherland particles on the circle 0„ are the eigenvalues of U, 
while the momenta p„ are the diagonal elements of Xi, specifically p„ = Bxn- The 
motion of the 0„ generated by the hamiltonian (^) is compatible with the evolution 
of the eigenvalues of U as it evolves in time according to (|4^). So the Sutherland 
model gives a holographic description of the quantum Hall state by monitoring the 
effective electron coordinates along X2, that is, the eigenvalues of U. 

The hamiltonians of the Sutherland and matrix model are equal and energy 
states map. The ground state is obtained by putting the particles at their static 
equilibrium positions. Because of their repulsion and the symmetry of the prob- 
lem, they will form a uniform lattice of points on the circle, as in pPj). 0o is the 
coordinate of the center of mass of the particles. Sound waves on this lattice cor- 
respond to small perturbations of the quantum Hall state. As the amplitude grows 
large, they become nonlinear nondispersive waves corresponding to holes forming in 
the quantum Hall particle distribution along the X-direction. In the limit they be- 
come solitons, representing isolated particles off the quantum Hall ground state ||3|. 
Further connections at the quantum level will be described in subsequent sections. 



5 Quantization of the cylindrical matrix Chern- 
Simons model 

5.1 Gauss law and quantization of the filling fraction 

We now proceed with the quantization of the above model. We use double brackets 
for quantum commutators and double kets for quantum states to distinguish them 
from matrix commutators and X-vectors. 
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The canonical structure of ( P5| ) implies the Poisson brackets 



{X^J,U,l} = -^SaU,, (51) 

This means that —BRX is the generator of right-rotations of the matrix U. Indeed, 
for any hermitian matrix e 

{-BRtT{eX),U} = iUe (52) 

Quantum mechanically X should also generate right-rotations of U. In the U- 
representation, then, it acquires the form 

As a result, TZ = —BRX satisfies the U{N) algebra: 

[[TZij,TZki]] = SalZkj - SkjTZii (54) 

with the obvious hermiticity condition Xjj = Xji. Similarly, the classical matrix 
C = BRUXU~^ generates left-rotations of U. It should therefore be ordered as 

A, = -t/..^ (55) 



and also satisfies the U{N) algebra (0). The sum of these operators Gjj = C + TI 
satisfies the SU{N) algebra {C and TZ have equal and opposite U{1) parts) and 
generates unitary conjugations of the matrix U. 

The components of are harmonic oscillators, satisfying 

[[^., = (56) 

The matrix = generates rotations of the vector \E', and must also satisfy 
the U{N) algebra (|^). It should therefore be ordered as 

(G*),, = (57) 

The Gauss law constraint then acquires the form 

{C + n + Gq,-9B)\phys)) = (58) 

The situation is similar to the planar case. Gu contains only symmetric products of 
the adjoint, with a number of boxes in their Young tableau (Z^v charge) a multiple 
of A^, kN {k integer). Gq, contains only totally symmetric representations with 
[/(I) charge equal to the number of boxes in the Young tableau. For the total 
representation to be in the singlet, as required by the Gauss law, Gu and G^r must 
be in conjugate representations and thus Gx^, must also have a number of boxes 
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kN . Moreover, the trace {U{1) charge) of Cj, must cancel NOB. So we obtain the 
quantization condition 

B6 = k , k = integer (59) 
Again, this is related to the level quantization of the noncommutative Chern-Simons 



action ^ and can also be attributed to a global gauge anomaly of the model 
4l[| . The above condition will lead to the quantization of the inverse filling fraction, 
as in Laughlin theory. We anticipate the actual result as z/ = + 1) = 1/n; the 
shift from /c to A; + 1 is a quantum correction. 

5.2 Quantum states 

The hamiltonian is 

It is the Laplacian on the group manifold U{N) (remember that X is essentially the 
momentum of U) , and is proportional to the common quadratic Casimir of C or TZ. 

Since the space U{N) is curved, we could add to the hamiltonian a term pro- 
portional to the curvature. Such terms can always arise through quantum ordering 
effects; a particular value, equal to 1/8 times the curvature, is singled out from 
conformal invariance. In our case the curvature is constant and, as we shall see, the 
addition of a constant term as above will make the spectrum especially simple and 
suggestive. 

Quantum states can be represented in terms of wavefunctions of U. A particular 
set of such wavefunctions diagonalizes H. Specifically, consider the matrix elements 
Raf}{U) of the matrix U in some irreducible representation of SU{N) R of dimension 
(Ir. By Schur's lemma, any arbitrary function of U can be expanded in terms of the 
above functions. 

Each matrix element Rap above is, in fact, an eigenstate of H. To see this, note 
that under arbitrary left- and right-rotations of U the states transform as 

Rap{V-^UW) = R-}^{V)R^s{U)Rsp{W) (61) 

So the multiplet RapiU), a, (3 = 1, ... transforms in the R representation under 
right-rotations of U and in the conjugate R representations under left-rotations of 
U . Since H is the quadratic Casimir of 7?. or £ we obtain 

HRap{U) = ^^C2,RRap{U) = EnRapiU) (62) 

So the spectrum of H consists of all quadratic Casimirs C2,_r, each with a degeneracy 
(P^ corresponding to the different matrix elements of Ra/siU). 

We still need to impose the Gauss law constraint. According to the discussion 
of the previous section, it stipulates that the states for U transform in a totally 
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symmetric representation Sk, with kN boxes, under conjugations of U. This means 
that we must pick the corresponding representation for Gu = C + IZ. Clearly the 
states RafiiU) transform in the Rx R representation under Gu- We must, therefore, 
decompose RxR into irreducible components and pick the symmetric representation 
Sk among the components. Each Sk corresponds to a unique physical state (the 
components within Sk are contracted in a unique way with the components of the 
\l/-representation Sk to give a singlet). This fixes the degeneracy of the eigenvalue 
Eji to the number of times Dji k that Sk appears in R x R. 

The above takes care of the SU{N) part of the wavefunction. We can always 
assign an arbitrary U{1) part by multiplying the wavefunction with (det Uy. If 
we want single-valuedness under the transformation U exp{i27r) U, which corre- 
sponds to rotations around the cylinder, we must have q = p/N with p an integer, 
corresponding to integer [/(l)-charge of the state under C or 71. Since the U{1) 
charge of R is the number of boxes, we can simply take the rows of its Young 
tableau to have either positive or negative length. 

We have reduced the physical spectrum of the model to pure group theory. In the 
present case, however, we can do better than that. Using standard Young tableau 
multiplication rules, it can be verified that in order to obtain Sk in the product 
RxR, the representation R must have lengths of Young tableau rows ij that satisfy 



ij+i 



>k, j = l,...N (63) 



For each such R, Sk is contained exactly once in R x R. 

The next interesting fact is that the quadratic Casimir of U {N) can be expressed 
in terms of the spectrum of free fermions on the circle |^ . Specifically, define the 
fermion 'momenta' 

+ 1 

p,=£, + ^— -J, j = l,...N (64) 
which satisfy pj > Pj+i- Then the expression for C2,r is 

N 

i=i 

where pj^ are the 'ground state' momenta, corresponding to the singlet representa- 
tion ij =0: 

Pj,o = — J (66) 

The final observation is that if we add to the hamiltonian a curvature term with 
coefficient 1/8, as mentioned in the beginning of the section, the 'ground state' 
term in C2,r exactly cancels. We are left, then, with the spectrum of free identical 
particles on the circle, but with an enhanced exclusion principle. Specifically 

^^ = E^. = ^Ep? (67) 
j j 
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where Ej are effective single-particle pseudo-energies, and the pj are single-particle 
pseudomomenta. Because of (0) and (|63D , the pj obey 

Pj - pj+i > k + 1 = n (68) 

We obtain the spectrum of free nonrelativistic particles of mass u/B on a circle of 
radius R, obeying exclusion statistics of order n = k + 1. This is the spectrum of 
Sutherland particles, which we have recovered entirely in the matrix model context. 

It is reasonable to interpret Ej as the quantum analogs of the eigenvalues of 
the potential {Buj/2)X^. This means that the positions of the electrons along the 
X-direction are 

1 , , 

^i = ^Pi (69) 

The ground state quasimomenta 

Pj,9s = ^^^-^^ — — C^O) 

form a 'Fermi sea' with distance n between successive momenta. The xj, then, 
correspond to evenly spaced electrons with a distance d = n/BR between them, 
and therefore an average density p = 1/{2ttR ■ d). The filling fraction then is 

'' = ^— (71) 

B n 

justifying the interpretation of n = + 1 as the quantized inverse filling fraction. 

= (the singlet sector) corresponds to free fermions, reproducing the fully filled 
n = 1 Landau level. 

Quasiparticle and quasihole states are identified in a way completely analogous 
to the planar case. A quasiparticle state is obtained by peeling a 'particle' from 
the surface of the sea (quasimomentum pi^gs) and putting it to a higher value pi > 
n{N — l)/2. This corresponds to an electron at position x ~ Pi/B along X in a 
state covariant under rotations of the cylinder. 

Quasiholes correspond to the minimal excitations of the ground state inside the 
quantum Hall tubular distribution. This is achieved by leaving all quasimomenta pj 
for j > r unchanged, for some integer r, and increasing all pj, j < r by one unit: 

N+1-2J . 
Pj = n , J >r (72) 

N+1-2J 

= n + 1 J <r (73) 

This increases the gap between pr and Pr+i to n + 1 and creates a minimal 'hole' at 
position X ~ Pr^gg/BR. As in the planar case, removal of one particle corresponds to 
the creation of n holes, and therefore the particle number of the hole is —q = —l/n = 



16 



—V. We again recover the quantization of the quasihole charge in fundamental units 
of 

qn^i^^- (74) 
n 

in accordance with Laughhn theory. 

Finally, we have center-of-mass {U{1)) excitations, achieved by shifting allpj by 
the same amount. This corresponds to translations of the electron state along the 
axis of the cylinder. 

6 Correspondence to Laughlin states 

The discussion in the last section demonstrates that there is a qualitative mapping 
between Laughlin states and matrix states. It is desirable to establish a more precise 
and explicit mapping between the two systems. 

The wavefunction of the electrons in the lowest Landau level is a function of 
two variables. Alternatively, we can use the reduced phase space representation, in 
which spatial coordinates do not commute and span a quantum phase space. The 
mapping between the two representation in the plane is through coherent states. 
Laughlin states can then be considered as particular restricted many-body Landau 
wavefunctions, or corresponding states in the reduced phase space. 

In order to establish the mapping between matrix and Laughlin states we will first 
define coherent states on a cylindrical phase space, then establish the correspondence 
with Laughlin states on the cylinder and finally map to matrix model states. 

6.1 Coherent states on a cylindrical phase space 

On a planar quantum phase space with coordinates X, Y satisfying [X, Y] = i/B we 
can define creation and annihilation operators a, as a = {X +iY)^jB/2. Coherent 
states, then, are defined as eigenstates of the annihilation operator a and represent 
states of minimum uncertainty. 

On a cylindrical phase space with one compact coordinate, say, Y = Y + 27ri?, 
the operator Y is multivalued and thus unphysical. The above creation-annihilation 
operators are therefore unphysical and we cannot use them to define coherent states. 
The operator exp[{X + iY)/ R], however, is single- valued and physical. We will then 
define coherent states \z) by the relation 

= e^/^\z) (75) 

We can easily find the expression for the wavefunction of \z). In the Y represen- 
tation, wavefunctions are periodic functions of Y and can be expanded in terms of 
the usual momentum states \n): 

(y|n) = e-^/«, n = -1,0,1,... (76) 
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In this representation X acts as {i/B)d/dY. A straightforward calculation shows 
that \z) is of the form 

— ^ rtz . . . . 

\z) ^ N^^e ^ \n) (77) 

n 

where is a normalization factor. Writing z = x + iy,we see that \z) is a state with 
X centered around x (although X has discrete eigenvalues) and Y centered around 
y (modulo 2TrR). 

It is convenient to choose the normalization N as 

TV = (27r3/2Bi/2^)-^/2g-^^'/^ = 7V„e-^^'/^ (78) 

Then it is straightforward to verify the completeness relation 

dzdz\z){z\ = l (79) 



Finally, the coherent state wavefunction of an n-state is 

• ny (n + BRx)^ 

{z\n) = NoC^ (80) 

or, defining w — exp{—x + iy) as the corresponding analytic coordinate on the 
cylinder, 

{w\n) = NoW''e~^''^~^ (81) 



6.2 Laughlin states on the cylinder 

We start by presenting the single-particle wavefunctions for the lowest Landau level 
on the cylinder. These are well-known, and can easily be found in the Landau gauge 
Ax = 0, Ay — —Bx. In a basis diagonalizing the momentum in the F-direction they 
become 

. ny B ( n \^ B 2 n'^ 

{x,y\n)j;^ = 'il)n{z^x + iy)^e^^~^y'^BR) = u'"e" 2"^ (82) 

They are 'stripe' states, exponential in the y-direction and gaussian in the x-direction 
with a center shifted by —n/BR. 

We immediately see the similarity with the coherent states in the reduced phase 
space. We conclude that Landau wavefunctions in terms of coordinates on the 
cylinder equal the corresponding coherent states on the reduced phase space. 

A'^-body Laughhn states on the cyhnder can be defined in a way analogous to the 
plane. There, the electrons were restricted to states containing the 'ground state' 
factor 

i^9s^\{{z,-zuTe-i^^\'^'^" (83) 

j<k 

For n = 1 this is the Slater determinant of the lowest angular momentum eigen- 
states, which reduces to the Vandermonde determinant for the variables Zi times the 
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(non-analytic) N-body oscillator ground state. For higher n the analytic (Vander- 
monde) part is raised to the n-th power, while the ground state part remains the 
same. 

A similar construction can be repeated on the cylinder, but with the momentum 
around the cylinder replacing the angular momentum. The corresponding Laughlin 
wavefunction would be 



The relevant single-particle states appearing in (R^) are, now, (p2|). The lowest 



i^9s = X{{w,-WkTe'-^^''^ (84) 



y-momentum appearing in (85) is 0, while the highes is n{N — 1), corresponding 
to positions in the a:-direction from to —n{N — 1)/BR. So this correspond to a 
tubular Hall state with area A ~ 2TTn{N — l)/B, corresponding to a filling fraction 
u = 1/n. 

To minimize the potential in x, the state should be centered around a; = and 
thus the above state should be shifted to the right by Ax = ^(A^ — 1)/2BR. This 
amounts to multiplying the wavefunction by the shift factor YljWj"'^'^ _ "YYi^ 
resulting, properly centered Laughlin wavefunction is 

= n(«.^)"e"^^^''^' , = te^ _ ' (85) 

For small distances the above wavefunction has the same behavior as the planar 
Laughlin state. It is interesting that its probability density also admits a plasma 
interpretation, as in the planar case. Indeed, the Coulomb potential on a torus takes 
the form 

V{x^ y) = - In (smh. — sinh — = :^ In ("^^ ~ ^"2") + 77 In {w^ — w~2~) (86) 



Therefore, the Vandermonde part of \'i)gs\ in (pSD will reproduce the exponential 
of the mutual Coulomb potential of particles on the torus, while the gaussian part 
represents the potential of a 'neutralizing' constant background charge distribution, 
in exact analogy to the planar case. 

The above ground state wavefunction can just as well be interpreted as the 
coherent wavefunction on the reduced phase space. Excited states above the ground 
state can be written by multiplying this wavefunction by any totally symmetric 
wavefunction of the Wi. Such wavefunctions can be uniquely expressed in terms of 
the Schur basis 

^c„c,...= (^E^'j l?"^")'" ^^^^ 

and the general state will be a product of Scx,c2... and tl}gs- 




V'ci.ca... = I E KM I E ^fc' 1 ■ ■ ■ n (^i - 



2 

e 2 -j 
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where we omitted the shift factor since, as all sjTiimetric functions, it can be repro- 
duced in terms of Schur functions. 



6.3 Mapping to matrix states 

The states of the matrix model can be explicitly written in a way analogous to the 
one for the plane ||2^. We will work in the U representation for the matrices and the 



oscillator representation for the \1'. Define a ground state wavefunction |0) which is 
the Fock vacuum of the oscillators and the singlet (constant) in U; that is, 

^,|0) =X,fc|0) = (89) 

Excited states can be obtained by applying \l'^'s and f/'s on |0). For the resulting 
state to be gauge invariant, all indices of Ujk and ^] must be contracted, either with 
each other or with the SU{N) antisymmetric tensor eJi-?2 - iJv_ 

The f/(l) gauge constraint (^), on the other hand, stipulates that each physical 
state should have exactly kN operators . The minimal way that we can contract 
the indices of the above \E'^'s is 

\gs) = (e^-^-^li^H^U),, ■ ■ ■ (vl/tf/^-i)^.^)' |0) (90) 

This can be considered the ground state of the system. The powers of U that appear 
must clearly be all different, and we chose them to span the values 0, 1, ... A^ — 1. We 
could have, instead, chosen the values —^^y^, ■ ■ ■ ^y^, which would have given the 
U part of the state a vanishing U{1) charge. This can always be done a posteriori 
by multiplying with det and we shall stay with \gs) above for simplicity. 

Other states can be obtained by multiplying with gauge invariant combinations 
of t/'s (no more are allowed). These are spanned by the Schur functions 

5ci,c2... = tr[/"Hrf/"2 . . . (91) 
and a complete basis for the matrix states is 

|ci, C2 . . . ) = tiWtiU''^ ■ ■ ■ {e^'-^''^^{^^U)j, ■ ■ ■ {^^U^~^)j^y |0) (92) 

The next step is to parametrize U in terms of diagonal and angular variables: 

U = V-^AV (93) 

with V an SU{N) matrix and A = diag{e'^^^} = diag{Wj} the eigenvalues of U. 
Each e-factor in the states (|92f) becomes 



e.i-.iv^;(vptt/)^.^...(v^tf/A^-i)^.^ 

{e''-'-V,,n ■ ■ ■ Vk^,^} {Wl ■ ■ ■ Wl-'} ■ ■ ■ i^^V~\,} (94) 
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Since the e tensor antisymmetrizes the indices jn, the indices /c„ appearing in Vk„j^ 
in the first bracket are also antisymmetrized and we obtain 

e^'-^^Vk,,,Vk,j, ■ ■ ■ Vk^j^ = e''^-'^ det V = e^-'^^ (95) 

(det V = 1 since V is special unitary). Defining 

X = V^ (96) 

the expression (|9^) becomes 

e'-'-W',^---W,l-'xl---xl (97) 



Since all fc„ are distinct, the product of the xl„ is simply Xixl' ' ' Xn- The remaining 
W^s with the e symbol reproduce the Vandermonde determinant. So the expression 
above becomes 

Ui^j-WJUx] (98) 

j<m j 

The Schur functions, on the other hand, become 

tiU^ = ^ (99) 
i 

Overall, the states of the matrix model (^) take the form 

C2 . . . ) = W^") f E ■■■UiW,- W^f [W x]\ |0) (100) 

The operators Xj defined above are also harmonic oscillators and they satisfy 
[[Xi'Xl]] = ^jk- So the oscillator state 

i^) = (^nx]j io) (101) 

appearing above has a norm independent of the matrix V which enters in the 

definition the Xy Since no other oscillator state can ever appear, the above state, 
used for calculations of matrix elements, is effectively independent of V and Wy 
The matrix V has, therefore, completely disappeared from the picture. 

When calculating matrix elements between states, we must integrate with the 
Haar measure \dU\ over the matrix V . This is 

\du\ = [dv] n i^i - ^A.r n (102) 

The integration [dV] over V produces a constant, since nothing depends on V. We 
are left with an integration over the eigenvalues (pj, with an additional Vandermonde 
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term coming from the measure. It is convenient to incorporate this measure into the 
definition of the states, so that matrix elements can be calculated with a fiat measure 
over the (pj. This introduces an additional power of the Vandermonde determinant 



in ( |100|) , shifting k to k + 1 = n. This is the origin of the renormalization of the 



filling fraction that we mentioned before. The final matrix states are, therefore, 

|ci,c2...)= [Y.wA lY.wA---I[iw,-w^r\n) (103) 

These are identical in form to the corresponding Laughlin states (^) upon mapping 
Wj = e*'^^ to Wj and identifying the ground state |0) with the bosonic gaussian factor 
in dH). 

The above provides a formal mapping between the states of the unitary Chern- 
Simons matrix model and Laughlin states on the cylinder, much like the one for the 



plane It should be stressed, however, that the above mapping is not unitary. 
Indeed, the wavefunctions ( |103| ) above are integrated with a fiat metric in (pj, while 
the Laughhn states ( pS] ) are integrated with the planar measure dxdy. Therefore, 
the norms and scalar products of the two sets of states are, in general, different. This 
mapping is, therefore, at best a qualitative one. The exact, unitary mapping between 
matrix model states and Laughlin states, planar or cylindrical, is still lacking. Some 
relevant results in the planar case are derived in . 



Note, further, that the above states ( |103| ) are not eigenstates of the eigenstates 
of the matrix model hamiltonian trX^. Eigenstates can always be constructed by 
forming appropriate linear combinations of states ( |103| ) with the same degree of 
homogeneity, it essentially amounts to constructing the characters of the appropriate 
allowed representations of U{N) (see the discussion in section 5.2). 



7 Outlook 

We have extended a previous proposal and presented a unitary matrix Chern-Simons 
model describing fractional quantum Hall states of electrons on the cylinder. 
The correspondence of the two systems was estabhshed as a map between states 
in each Hilbert space. As in the planar case, the quantization of the inverse filling 
fraction and of the quasihole charge are straightforward consequences of the quantum 
mechanics of this model. We also stressed that the classical value of the inverse filling 
fraction is shifted quantum mechanically by one unit. This can be equivalently 
viewed as a renormalization of the Chern-Simons coefficient, as a group-theoretic 
effect or as a result of the nontrivial measure of the model. 

We further pointed out that this model, and therefore also the two-dimensional 
quantum Hall system, is described holographically in terms of a one-dimensional 
system, the so-called Sutherland integrable model of particles on the circle. 

The correspondence between states and operators of the matrix and quantum 
Hall systems is still an open issue. In principle, such a correspondence is guaranteed 
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to exist, since the two Hilbert spaces have the same dimensionahty, as demonstrated 
in the last section. For it to be useful, however, it should be such that operators in the 
quantum Hall system map into explicit, simple matrix model operators. The hope is 
that the operators in the matrix model language will not involve explicitly the filling 
fraction (unlike, e.g., the hole creation operators in the second-quantized fractional 
quantum Hall system), and thus will describe the properties of these systems in a 
more universal way. This would also open the road for the calculation of relevant 
quantities, such as correlation functions, in the matrix model formulation. This is 
a most important issue. 

There are clearly many other open questions. Incorporating the spin of the 
electrons and identifying skyrmion-like configurations is an obvious next step. Most 
intriguing, however, is the question of a possible phase transition of the quantum 
Hall system at small filling fractions. Numerical simulations suggest that at z/~^ ~ 67 
Laughlin electrons form a Wigner crystal instead of an incompressible fiuid. This is 
based on the properties of the Laughlin wavefunctions and does not seem to hinge 
on the specific dynamics of the electrons beyond what is already encoded in the 
wavefunctions themselves. The corresponding one-dimensional Sutherland system 
does not exhibit any such phase transition. This does not guarantee, however, that 
two-dimensional quantities calculated in the context of this model would not exhibit 
nonanalytic (or at least crossover) behavior in the filling fraction, signaling a phase 
transition. This intriguing possibility is the subject of further research. 
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